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Phase transitions in isotropic quantum antiferromagnets are described by an 0(3) nonlinear quan¬ 
tum field theory. In three dimensions, the fundamental property of this theory is logarithmic scaling 
of the coupling constant. At the quantum critical point the coupling asymptotically vanishes and 
the quasiparticles become free. This logarithmic decay of the coupling constant has never been ob¬ 
served. In this paper, we derive finite temperature properties of the field theory and use our results 
to analyze the existing data on the real antiferromagnet TlCuCla. Including finite temperatures 
into the theory, we find that agreement between theory and experiment is sufficiently sensitive to 
unambiguously identify the asymptotic decay of the coupling constant. We also comment on the 
unique possibility to study Landau pole-physics in quantum magnets. 
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Asymptotic freedom plays a crucial role in Quantum 
Chromodynamics. The freedom, which is due to non- 
abelian gauge fields, means a logarithmic decay of the 
coupling constant at high energies. Ultimately at infi¬ 
nite energies particles do not interact, this is ultraviolet 
asymptotic freedom mm- In three dimensional (3D) 
non-gauge quantum field theories as well as in abelian 
gauge theories, the coupling constant decays logarithmi¬ 
cally at low energies [3]. However, usually this decay is 
terminated because of a low energy cutoff. For exam¬ 
ple in Quantum Electrodynamics the cutoff is due to the 
rest energy of the electron. The low energy logarithmic 
decay can be observed only at a quantum critical point 
(QCP) where the cutoff energy is zero. It is known to the 
theory of quantum magnetism that certain antiferromag¬ 
netic systems can be driven to a QCP, and that at such a 
point the characteristic low energy magnetic excitations 
become non-interacting. At the QCP they experience 
“infrared asymptotic freedom”. Quantum antiferromag¬ 
nets in the vicinity of a QCP therefore provide a perfect 
testing ground for such behaviour. 

The low energy logarithmic decay of the coupling con¬ 
stant at a QCP can, in principle, be observed at zero tem¬ 
perature. The zero temperature case is well understood 
theoretically, a 3D quantum antiferromagnet at zero tem¬ 
perature is equivalent to a 4D classical antiferromagnet 
at finite temperature. The 3D QCP corresponds to the 
Neel transition in the 4D case. Thermodynamic quanti¬ 
ties scale as powers of the running coupling constant [4]. 
Unfortunately, the existing zero temperature experimen¬ 
tal data are not sufficient to pin down the logarithmic 
scaling. On the other hand combined zero and nonzero 
temperature data on TlCuCls provide an excellent 
opportunity to search for fingerprints of asymptotic free¬ 
dom at the QCP. To perform this search we develop a 
theory of the QCP which accounts for both quantum and 
thermal fluctuations. After that we compare the theoret¬ 
ical predictions with experimental data. The comparison 
unambiguously indicates a logarithmically running cou¬ 



FIG. 1: The phase/energy diagram of TlCuCls The ver¬ 
tical panel shows the pres sure-temperature phase diagram, the 
Neel temperature curve separates magnetically ordered and 
magnetically disordered phases. The light red band around 
the Neel curve indicates the region of dimensional crossover. 
The horizontal panel shows both the triplon gap At in the 
paramagnetic phase and the Higgs magnon gap Ah in the an¬ 
tiferromagnetic phase versus pressure at zero temperature. 


pling constant. 

The phase diagram of the dimerized 3D quantum an¬ 
tiferromagnet TlCuCls is shown in the vertical panel of 
Fig.m The disordered quantum state consists of an ar¬ 
ray of spin dimers (spin singlets), and the ordered quan¬ 
tum state has a long range Neel order as illustrated in 
Fig. The Neel temperature curve (red line) separates 
ordered and disordered phases, the QCP is indicated by 
the yellow dot. 

Excitations in the disordered phase, triplons, are 
gapped. These are triplet excitations of spin dimers, see 
Fig. IJi. There are two different kinds of excitations 
in the ordered phase, gapless Goldstone excitations and 
the gapped longitudinal Higgs excitation, they are illus¬ 
trated in Eig. [^ & c. The horizontal panel in Eig. 
displays excitation gaps versus pressure at zero temper¬ 
ature. Overall the experimental data m\ provide the 
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FIG. 2: Excitations of a dimerized quantum antiferromagnet. 
Panel a illustrates the triple degenerate gapped triplon excita¬ 
tions. Panel b illustrates the gapped longitudinal (Higgs) ex¬ 
citation. Panel c illustrates the quantum phase transition; the 
strength of the interactions in either phase is depicted by the 
steepness of the ‘well’. Within the ordered phase, the ‘Mexi¬ 
can hat’ potential has a flat direction which supports the gap¬ 
less Goldstone excitations (red arrows). Precisely at the QCP 
(dashed line), all directions flatten - the Higgs and triplon 
excitations become gapless and non-interacting i.e. asymptot¬ 
ically free. 


following information, (i) The Neel temperature versus 
pressure, (ii) The magnetic excitation gap in the disor¬ 
dered phase for various temperatures and pressures, (iii) 
The Higgs magnon excitation gap in the antiferromag¬ 
netic phase for various temperatures and pressures, (iv) 
The magnetic excitation width (lifetime) for various tem¬ 
peratures and pressures. In our analysis we do not use 
fully the width data, the width is used only to indicate 
the dimensional crossover region around the Neel tem¬ 
perature. However, we fully use the data from points 
(i), (ii) & (iii). In the analysis we neglect small spin- 
orbit anisotropy which leads to a small gap in one of the 
“Goldstone” modes in the antiferromagnetic phase. 

The quantum phase transition (QPT) between ordered 
and disordered phases is described by the effective field 
theory with the following Lagrangian iHHn], 

^=^0^,08^^0-^17110^-^ao[0^f. ( 1 ) 

The vector field (fl describes the staggered magnetisation, 
the index fi enumerates time and three coordinates. The 
QPT results from tuning the mass term, mg, for which 
we take the linear expansion mg (p) = 7 ^ {pc — p ), where 
7 ^ > 0 is a coefficient and p is the applied pressure. 
Varying the pressure leads to two distinct phases; (i) for 
p < Pc we have rug > 0 , and the classical expectation 
value of the field is zero = 0. This describes the 
magnetically disordered phase, the system has a global 
rotational symmetry, and the excitations are gapped and 
triply degenerate. These excitations are referred to as 
“triplons”. (ii) For pressures p > Pc we have mg < 0, 
and the field obtains a non-zero classical expectation 
value ip‘1 = This describes the magnetically or¬ 


dered, antiferromagnetic phase. Varying mg from posi¬ 
tive to negative spontaneously breaks the 0(3) symmetry 
of the system. In the broken phase there are two gapless 
transverse (Goldstone [ 12 ]) excitations, and one gapped 
longitudinal (Higgs) excitation. One easily recovers the 
known relation; Higgs gap/triplon gap= \/2 [13], explic- 
itly At{p) = mo{p) and A^(p) = V 2 |mo(p)|. 



FIG. 3: Diagrams for the vertex F and self-energy T. 


The above analysis does not account for quantum or 
thermal fluctuations. All fluctuations considered in the 
present paper originate from the vertex and self-energy 
diagrams shown in Fig. The vertex corrections result 
in the running coupling constant a a, see e.g. Ref. [4] or 
the Supplemental Material, 

^ _ _ ^0 _ 

A l + ^ln(Ao/A)- 

Here A is the energy/momentum scale, and Ag is the nor¬ 
malization point, q^Aq = Qfg. Eq. (H has been obtained 
within the single loop renormalization group which im¬ 
plies that ao/Sn 1. At the same time the logarith¬ 
mically enhanced denominator in § can be arbitrarily 
large. Note that the normalization point Ag can be ar¬ 
bitrary, generally it is not equal to the ultraviolet cutoff 
related to the lattice spacing. Eq. has a pole at 
A = Ap = Age^^ /iiao^ This is the famous Landau pole 
much debated in quantum field theory [3]. Physics in the 
vicinity of the Landau pole remains a controversial issue 
in spite of 60 years of theoretical studies. Remarkably 
quantum magnets can shed light on the problem, we will 
return to this point later. The running coupling constant 
at zero temperature versus pressure is plotted in Eig. 
This curve is extracted from available experimental data. 
The coupling constant vanishes at the QCP, indicating 
the asymptotic freedom. Eig. represents one of our 
central results, below we explain how we obtain it. 

At zero temperature, equations for the running mass 
and for the running staggered magnetization are well 
known [4], (see also the Supplemental Material) 


rn^iP,^) = l'^{Pc-p) 
l^{p-Pc) 


p0p,h) = 




ao _ 
<Ao 

(Aa 


( 3 ) 

( 4 ) 


To find actual values of the gap in the disordered phase 
one has to solve Eq. © at p < Pc with At = A = m. 
To find the Higgs gap in the ordered phase one has to 
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Pressure [kbar] 

FIG. 4: Zero temperature running eoupling eonstant versus 
pressure in TlCuCls. The eonstant vanishes at the QCP (yel¬ 
low point). 


solve Eq. ^ at p > Pc with Ah = A = ^/2\m\. The 
relation Ah /At = a/ 2 [13] remains valid with logarithmic 
accuracy. 

We need to extend the theory to nonzero temperatures. 
Our goal is to find excitation spectra, therefore we can¬ 
not use the imaginary time Matsubara technique. We 
need to work with real frequencies at a nonzero tempera¬ 
ture. Generally, there is not a regular diagrammatic tech¬ 
nique which allows to calculate real frequency Green’s 
functions at nonzero temperature. Fortunately, in the 
present case such a calculation is possible within standard 
techniques. This possibility comes from the two follow¬ 
ing observations, (i) Multi-loop logarithmic corrections 
are universal, they are not sensitive as to whether fre¬ 
quency/energy is real or imaginary, (ii) The leading in a 
correction which contains powers of temperature comes 
only from the self energy diagram shown in Fig. Cal¬ 
culation of this diagram does not cause problems since 
while it depends on temperature, it is frequency indepen¬ 
dent. Still, there is a minor complication related to point 
(ii). The complication is due to temperature broaden¬ 
ing (quasipraticle lifetime) because of scattering from the 
heat bath of magnons. Below we explain how we address 
this issue. To be specific let us consider the triplon gap 
in the disordered phase. Calculation of the self-energy F, 
Fig. i gives the following answer (see the Supplemental 
Material) 


a|(p,t,a) 


7^(Pc-p) 


ao 


11 _^ 

k 


i/wfe 


e T 


- 1 


• (5) 


At zero temperature the second term in the right hand 
side is zero and Eq. ^ becomes identical to Eq. 

The value of the lower logarithmic cutoff in (§ is obvious, 
A = max{At,T}. The triplon dispersion is harder. The 
naive formula uuk = y/k‘^ + is incomplete because at 
small k and close to the Neel temperature where At 
0 the line width Tt (temperature broadening) becomes 


larger than the gap. Physically the inequality Ft > At 
is an indication of the dimensional crossover, AD 3D. 
Sufficiently close to the Neel temperature, critical indices 
take the 3D classical values. To fix the problem we take 
ujk = + r|, this is a standard way to describe 

a damped harmonic oscillator, see e.g. Ref. [5|. Of 
course the modified dispersion is not sufficient to fully 
describe the dimensional crossover, but it is sufficient for 
the purposes of the present work. The line broadening we 
take directly from experiment, T^ = ^T, where ^ « 0.15 
0. Solution of Eq. ^ with At = 0 gives the Neel 
temperature as function of pressure, Th{p)- 

One can also approach the Neel temperature from the 
ordered phase, see Supplemental Material. In this case 
there are two Goldstone and one Higgs mode. Equation 
for the Higgs gap is similar to 


A|f(p,r,A) =2<7^(p-pc) 


(Aa 


ao 


(6) 


^ eT -1 “ 


^ 1 

e ^ — 1 


Again, A = max{Ai:/,T}, uuk = 

Th = (T. The Neel temperature determined from the 
condition At = 0, Eq. ([^, must be identical to that de¬ 
termined from Ah = 0, Eq. From here we conclude 
that broadening of the Higgs mode is larger than that of 
the triplon, C 0.3 compared to ^ ~ 0.15. The larger 
broadening is consistent with the data [5] . We note that 
the critical exponent of the magnetisation in Eq. 
is identical to that found for the Neel temperature by 
solving Eq. or (H- This agrees with the latest Quan¬ 
tum Monte Carlo simulations on the three dimensional 
dimerised antiferromagnet m- 

Now we are fully armed to perform fits to the experi¬ 
mental data. As the normalisation point we take Aq = 1 
meV. We remind the reader that this choice is arbitrary, 
one can always use a different normalisation point with 
an appropriate rescaling of the coupling (Aq. There are 
three fitting parameters, the critical pressure Pc, the co¬ 
efficient 7 ^ in the pressure dependence of the bare mass, 
and the coupling constant ao- Points in panel a of Fig. 
show experimental values of the triplon and Higgs gaps 
for various pressures at T = 1.85 K. Points in panel b 
of Fig. show Neel temperatures for various pressures. 
The data are taken from Refs mi- 

Solid curves in both panels show fit of the data with 
Eq.’s (H, (§ and <©• Values of the fitting parameters 
are 

Pc = 1.01 kbar, 7 = 0.68 meV/kbar^/^, ^ = 0.23 . (7) 

Sir 

It is worth noting that while T = 1.85K is a pretty low 
temperature, the temperature corrections in Eq.’s (§ and 
§) are not negligible. In the present work we set the 
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FIG. 5: Panel a: triplon and Higgs gaps versus presssure at 
temperature T=1.85K. Points show experimental data from 
RefE- Panel b: Neel temperature versus pressure. Points 
show experimental data from Ref m- In both plots the solid 
and dashed eurves are quantum field theory fits with and with¬ 
out aeeount of the running eoupling eonstant, respeetively. 

triplon speed equal to unity. If one restores three dif¬ 
ferent speeds along three different principal directions 
of the lattice, ci, C 2 , C 3 , then Eq. 0 is changed to 
ao/( 87 rciC 2 C 3 ) = 0.23. This value is close to the value 
ao/( 87 rciC 2 C 3 ) = 0.21 obtained in Ref. [10] from the 
lifetime of the Higgs mode. The work m did not ac¬ 
count for the running coupling constant, however ex¬ 
perimentally the major contribution to the lifetime data 
comes from Higgs magnons with energy of about ImeV. 
This energy is taken as the normalisation point in the 
present work. Coincidence of the energy scales explains 
the very close agreement between the accurate result of 
the present work and the approximate result of Ref. [T0| . 

We stress that the coupling constant significantly 
changes along the fitting curves in Fig. To illus¬ 
trate this change we present Fig. which shows how 
the constant runs along the Neel temperature curve. A 
similar run is shown in Fig. where the coupling con¬ 
stant is plotted versus pressure at zero temperature. In 
this case the infrared cutoff in Eq. 0 for q^a is equal to 
triplon/Higgs gap at zero temperature. Position of the 
Landau pole follows from the known value of the cou¬ 
pling constant, Ap = /n^o ~ 3.5meV. This energy 


FIG. 6: Running eoupling eonstant versus pressure along the 
Neel temperature eurve. Unlike Fig. 4 where the temperature 
is zero, in this ease T = Tn{p)- Aeeording to Eq.’s ^ and 
0. the infrared eutoff in Eq. ^ for oa is A = Tn{]^- The 
QCP is again marked by the yellow dot. 


is higher than the experimentally studied regime and is 
comparable with expected ultraviolet cutoff related to 
the dispersion along the third axis, see discussion in Ref. 
HO]. Experimental and theoretical studies in this energy 
range can shed light on the Landau pole physics and on 
the expected dimensional crossover in TICUCI 3 . Alter¬ 
natively the Landau pole physics can be addressed in 
Quantum Monte Carlo studies of dimerized spin-lattice 
models [HIISI. 

Our central goal is to pin down the running coupling 
constant which vanishes at the QCP giving rise to asymp¬ 
totic freedom. To check this statement we have also per¬ 
formed a fit of data with fixed coupling constant. For this 
purpose we use the same Eq’s. 0 and 0 where we set 
= 0 ^ 0 - The best fitting parameters are Pc = l.Olkbar, 
^ = 0.64meV/kbar^/^, ao/ 87 r = 0.16. Corresponding 
fitting curves are shown in Fig. [^by dashed lines. It is 
seen from the quality of the fits that the running coupling 
plays a crucial role in describing the static properties of 
the system; the analysis clearly demonstrates the running 
coupling constant. 

Asymptotic freedom is a prominent physical phe¬ 
nomenon. It is the remarkable experimental control of 
the quantum antiferromagnet TICUCI 3 that has allowed 
the present work to identify for the first time the loga¬ 
rithmic decay of the coupling constant. More generally, 
with such remarkable experimental control, TICUCI 3 and 
other quantum antiferromagnets provide an ideal play¬ 
ground for studies of the Landau pole physics and many 
other nontrivial quantum phenomena. 

We thank Bruce Normand, Anders Sandvik and 
Yaroslav Kharkov for many useful discussions. We also 
thank Christian Riiegg for his comments regarding the 
decay widths. 
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A. Running Coupling Constant 


write the two point function as 


The four point vertex in Fig. (3) is calculated to second 
order in a (with a Euclidean metric) 


rW = a-lla- 


L 


d^k 1 
A (27r)4 


lla^ f Ac 

= ct- In -r 

V a 


(A.l) 


The infrared cut-off, A, is given by the mass gap, or the 
temperature scale. We use a Callan-Symanzik equation 
to find the Beta function 


d R( \ ^ 

d\n{Ac/A) ^ da 


r(4) = 0 


P{a) 


da 


(iln(Ao/A) 

q^a = 


IW 

87r^ 

lla^ 

87r^ 


ao 


^ln(Ao/A) 

(A.2) 


p(2) _ , 


5(aA 


L 


d^k 


A (27r)3 2v^P~+~^ 


m 


5q!ato2 (Ac 


(BA) 


We use the Callan-Symanzik equation to find the (mass) 
Beta function 


0 = 


d 


_dln(Ae/A) 


+ /^m(A) 


d 

dm? 


=> /3m (A) = 

87r2 

dm? 

ba^m? 

(iln(Ao/A) 

87r2 


(—] 



lllj 

l + ^ln(Ao/A) 

d\n{m^) 

(-) 

, ^“0 

(iln(Ao/A) 

(llj 

l + ^ln(Ao/A) 


2 2 / 

= ^0 — 
\<Ao 


p(2) 


(B.4) 


where Ac is some momentum cut-off such as the inverse 
lattice spacing, while Aq is the ‘normalization’ scale or 
point. 

B. Self-Energy in the Disordered Phase 


Including non-zero temperatures does not change the 
form of the running coupling nor mass eq.’s ( |A.2[ |B.4[ ), 
but it does shift the infrared cutoff from m{p) 
A=Max{A,(p,T),T} 

After accounting for how the coupling terms and A 
depend on the scale, discussion given below, we find that 
the gap takes the form 


Approaching from the disordered phase, the first pertur¬ 
bative correction to the triplon gap comes from the one 
loop self energy 


E(A,T) =5aA V — 

k 


e T 


- 1 


(B.l) 


= 5q;a 


/ 


d^k 1 


+ 5q;a 


/ 


d^k 1 


(27r)3 Wfc (e^ - 1)’ 


The coupling constant coefficient is the running coupling 
a A, since the two point corrections are multiplicative 
with the four point vertices. With these corrections the 
triplon gap becomes dependent on both p and T 


A2(p,T)=m2(p) + E(A,r). (B.2) 


^tiP,T,A) 




p) 


<AA 

<A 0 


_ 5 _ 

11 


+ 5(aA ^ ^ 

k 


1 1 

— 1 
(B.6) 


C. Self-Energy in the Ordered Phase 


The ordered phase is induced by the spontaneous break¬ 
down of the 0(3) symmetry when p > Pc, 8iS discussed 
in the introduction. By analogy with the result for the 
mass gap in the disordered phase, we find that the Higgs 
mass in the ordered phase takes the form 


Al{p,T,A)=2 


R{p-Pc) 


aA 


_ 5 _ 

11 


(C.l) 


The first term in the self energy eq. renormalizes the 
bare mass term ttIq, such that mQ-\-baA f 2 ^ ^A 
has logarithmic dependence on the energy scale A. After 
RG, this part results in Eq.(3). The second term, or the 
‘temperature perturbation’, only contributes to the log¬ 
arithmic running via its influence on the infrared cutoff. 
To make these statements more clear, consider zero tem¬ 
perature such that only the first term contributes. We 


- 2 q(a 

k 


1/k 

(e^ - 1 ) 


3<aa 

k 


1 /^fc [ 


Again, A = max{Aff,T}, and the scale dependence is 
of course the same as in the disordered phase. It is a 
delicate task to calculate the self energy contributions to 
the Higgs gap, since within the ordered phase our calcu¬ 
lations at each order in A must preserve the Goldstone 
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theorem. The Goldstone theorem is a direct result of 
the remaining 0(2) symmetry. We outline the proce¬ 
dure here; In the Lagrangian, the field (p = (tt, (pc + cr) 
is shifted such that the minimum of the potential is pci 
and the field oscillations about this shifted minimum are 
the two Goldstone modes tt and the gapped Higgs mode 


We can write an effective potential, V, from the non¬ 
derivative terms of the Lagrangian expanded about the 
the minimum pc 

V = -i|m^|(^,99c +cr)^ + + (C.2) 

The following two conditions must simultaneously hold 
true to ensure that pc is indeed the minimum of the po¬ 
tential, and that to any order in A, the perturbations re¬ 
spect the 0(2) symmetry and so preserve the Goldstone 
theorem 


dV 

dp 


= 0 , 


and 


dTT^ 


= 0 . 


(C.3) 


Since we have already obtained the universal scale de¬ 
pendence of a A and tua, we do not need to repeat the 
Gallan-Symanzik, RG procedure. We just outline how 
the thermal perturbations are to be treated. Gomputing 
the thermal loops explicitly we obtain the first expression 


dV 

dp 


= C^APc 




E l 

k 


1/k 


1 ) 


■ 3aA 2^ — 


(e' 


= 


RaI 

q;a 


1 ) 

1/k 


= 0 


(C.4) 


-2E7^-3E74 


l/uJk 


{eT - 1 ) 


(e“ 


- 1 ) 

(G.5) 


where the thermal corrections are split into two separate 
contributions. This is because one type comes from the 
one loop self-energy with a Higgs propagator, and the 
other with a Goldstone propagator. The first summation 
accounts for loops with massless Goldstone propagators, 
while the second accounts for loops with massive Higgs 
propagators, so that Cc;| = + A//(p,T)^. We can now 

find the Higgs gap using the result from Eq.(G.5). Di¬ 
rectly computing the one loop corrections to the Higgs 
gap, we find 


A% = {SaApc - kil 


k (eT - 1) (e T - 1) J 


= 2|mA|"-4aAE74^-6«AE7^ 

k (eT - 1) k (e - 1) 

(C.6) 


1 /wfc 


and we have used Eq.(G.5) in passing from the first to 
second lines. We see that A‘jj = 2aAp\ + 0{a‘^) which is 
equivalent to the result for the Higgs gap Eq.(6). 


D. Neel Temperature 


Approaching from the disordered phase, we calculate 
the Neel temperature by solving Eq.(5) for At(p, T/v) = 0 


TNipf 


A{p-Pc) 

5«oEy(5fe) 


0^0 

aA 


11 


(D.l) 


where ujy = -1 - (T/Ta)^ = The fit T = 

was discussed in the main text. 

Similarly, we can approach from the ordered phase 
and calculate the Neel temperature by solving Eq.(6) for 
Ah{p,Tn) = 0 


Tn{p)^ = 


l‘^{p-Pc) 




1 /y 
y {eV-l) 


0^0 

aA 


(D.2) 


here Cjy = f^e two terms in the denom¬ 

inator are due to the Higgs and Goldstone self-energies. 


Since the phase transition is of second order, Eq.’s (D.l) 


and (D.2) are equivalent. Glearly without account of the 


finite line width T, the equations are identical. As dis¬ 
cussed in the main text, we approximate T^ = and 
Th = from the experimental data. This approxima¬ 
tion only becomes important in the vicinity of the phase 


transition Equating Eq.’s (D.l) and (D.2), 

we find how ^ and ( are related at the phase transition. 
Again, we take ^ ^ 0.15 and ( ^ 0.3. 
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